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ABSTRACT 

Recently. iFossati et al.l observed that the UV transit of WASP-12b showed an early ingress compared 
to the optical transit. We suggest that the resulting early ingress is caused by a bow shock ahead 
of the planetary orbital motion. In this Letter we investigate the conditions that might lead to the 
formation of such a bow shock. We consider two scenarios: (1) the stellar magnetic field is strong 
enough to confine the hot coronal plasma out to the planetary orbit and (2) the stellar magnetic 
field is unable to confine the plasma, which escapes in a wind. In both shock capable of 

compressing plasma to the observed densities will form around the planet for plasma temperatures 
T < (4 — 5) x 10 6 K. In the confined case, the shock always forms directly ahead of the planet, but in 
the wind case the shock orientation depends on the wind speed and hence on the plasma temperature. 
For higher wind temperatures, the shock forms closer to the line of centers between the planet and 
the star. We conclude that shock formation leading to an observable early UV ingress is likely to be a 
common feature of transiting systems and may prove to be a useful tool in setting limits on planetary 
magnetic field strengths B p . In the case of WASP-12b, we derive an upper limit of about B p = 24 G. 

Subject headings: planet-star interactions — planets and satellites: individual (WASP-12b) — planets 
and satellites: magnetic fields — stars: coronae — stars: individual (WASP-12) — 
stars: winds, outflows 



1. INTRODUCTION 

WASP- 12b is among the largest transiting planets 
discovered so far. First identified in an o ptical pho- 
tometric transit survey (jHebb et al.l I2009L hereafter, 
H09), its mass is M p = 1.41 Mj and radius R p = 
1.79 Rj, where Mj and Rj are the mass and radius 
of Jupiter, respectively. After its discove ry, additional 
observ at ions have been acquired ( e.g.. IFossati et all 
2010aHbl: lLopez-Morales et all 120101: iHusnoo et alJl201fi 
Campo et al.l 12010). Interestingly, transit observations 
in the near-U V revealed a longer t ransit duration than 
in the optical (|Fossati et al.ll2010al ). While the time of 
the egress occurs almost simultaneously for both set of 
observations, the ingress of the transit is first seen in the 
near-UV wavelength range. This asymmetric behavior 
has been explained by the presence of asymmetries in 
the planetary atmosphere. 

Close-in giant gas planets are rather in- 
flated and most have developed a n ex os phere 
(HP 209458b: IVidal-Madiar et all 120031 120081 : 
Ehrenreich etaLl 120081: iLinskv et al.l l2010l: HD 189733b: 
Lccavclier Dcs Etangs et al. I2010L CoRoT-Exo-lb : 
Barge et al.l 120081: WASP^12b: IFossati et al.l l20T0al) 
that can fill or even overflow the planet's Roch e lobe 
(IGu et al.ll200l ILi et al.l[2010t llbgui eT^TI [2010^1 . This 
may result in mass transfer through a Lagrangian 
point to the star that could cause an asymmetry in 
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the appearance of the t ransiting planet-star system as 
seen from the Earth ()Lai. Helling, fc van den Heuvell 
120101 hereafter, LHV10). Asymmetri es could also be 
produ ced by cometary tails. However, lEhrenreich et all 
(|2008ft demonstrated for HD 209458b that a radiation- 
driven cometary tail would produce a late egress of the 
planeta ry transit light curve, instead of an early ingress. 
ILHVIOI investigated the formation of a bow shock 
around the planet due to the interaction of the planet's 
magnetosphere with a stellar wind as the cause of the 
early ingress. Assuming a typical solar wind m ass-loss 
rate and adopting solar wind properties, ILHV10I derived 
the wind velocity at the dist ance of WASP - 12b using a 
thermally driven wind model ((Parker 1958J. They found 
that the wind is still subsonic at WASP- 12b orbital 
distance, concluding that no bow shock is expected 
to form in the planet- wind interaction zone. Here, we 
demonstrate that a bow shock can actually be formed 
around the planet if the relative azimuthal velocity 
between the planetary orbital motion and the ambient 
medium is taken into account. 

WASP-12b orbits its host star (a late-F main-sequence 
star, with mass Af* = 1.35 M Q and radius i?» = 1.57 Rq) 
at an orbital radius of i? rb = 0.023 AU — 3.15 with 
an orbital period of P orb = 1.09 d (|H09t ). Due to its close 
proximity to the star, the flux of coronal particles im- 
pacting on the planet comes mainly from the azimuthal 
direction, as the planet moves at a Keplerian orbital ve- 
locity of u K = (GA/»/i? orb ) 1/2 ~ 230 km s" 1 around the 
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star. Therefore, stellar coronal material is compressed 
ahead of the planetary orbital motion, possibly forming 
a bow shock ahead of the planet. If such compressed ma- 
terial is optically thin, the planetary transit light curve 
is symmetrical with respect to the ingress and the egress 
of the planet. Indeed, this is the case w hen the transit is 
observed at optical wavelengths (H09). However, if the 
shocked material ahead of the planet can absorb enough 
stellar radiation, the observer will note an early ingress 
of the planet in the stellar disk, but no difference will 
be seen at the time of egress, as the shocked material is 
present only ahead of the planetary motion. 

Here, we investigate under which conditions the inter- 
action of a planet with the stellar coronal plasma could 
lead to the formation of a bow shock ahead of the plan- 
etary orbital motion, and therefore, explain the early 
ingress observable in the near-UV. Our shock model is 
described in Section [2] Although we know the orbital 
radius of the planet, we do not know if at this radius the 
stellar magnetic field is still capable of confining the hot 
gas of the stellar corona, or if this gas is escaping in a 
wind. Therefore, we investigate the validity of our model 
with respect to different stellar coronal conditions in Sec- 
tions [3] and 0] Conclusions are presented in Section [5] 
2. THE SHOCK MODEL 

A bow shock around a planet is formed when the 
relative motion between the planet and the stellar 
corona/wind is supersonic. The shock configuration de- 
pends on the direction of the flux of particles that ar- 
rives at the planet. We illustrate two different limits 
of the shock configuration and an intermediate case in 
Figure [IJ where 9 is the deflection angle between the az- 
imuthal direction of the planetary motion and n, a vector 
that defines the outward direction of the shock. As seen 
from the planet, — n is the velocity of the impacting ma- 
terial. 

The first shock limit, a "dayside-shock" , occurs when 
the dominant flux of particles impacting on the planet 
arises from the (radial) wind of its host star. For in- 
stance, the impact of the supersonic solar wind forms a 
bow shock at the dayside of Earth's magnetosphere (i.e., 
at the side that faces the Sun). This condition is illus- 
trated in Figure QJa) and is met when u r > c s , where 
u r and c s are the local radial stellar wind velocity and 
sound speed, respectively. 

A second shock limit, an " ahead- sho ck" , occurs when 
the dominant flux of particles impacting on the planet 
arises from the relative azimuthal velocity between the 
planetary orbital motion and the ambient plasma. This 
condition is especially important when the planet orbits 
at a close distance to the star, and therefore, possesses a 
high uk- In this case, the velocity of the particles that 
the planet 'sees' is supersonic if Am = \uk — u v \ > c s , 
where u v is the azimuthal velocity of the stellar corona. 
This condition is illustrated in Figure [ljb). 

For intermediate cases, both the wind and the az- 
imuthal relative velocities will contribute to the forma- 
tion of a shock around the planet (Figure QJc)) and the 
deflection angle 9 is given by 



= atan h r . (1) 

\\u K - u v \J 

In general, for a planet orbiting its host star at a close dis- 



tance, the stellar wind is still accelerating and subsonic 
[u r < c s ). In this case, conditions for an ahead-shock 
will more probably be met. 

Here, we adopt the estimates of ILHVIOI to derive the 
minimum density of the stellar plasma at the orbital ra- 
di us of WASP- 1 2b. F r om the observations presented 
in iFossati et al.l ()2010aD . ILHVIOI estimated the column 
density of the absorbing gas around the planet to be 
> 1.4 x 10 13 cm~ 2 . For that, they assumed r = 1 in the 
absorption lines of Mg. Because the time resolution of 
the transit measurement is sparse, it is difficult to place a 
firm constrain on the thickness of the absorbing material, 
which could either extend all the way to the planetary 
surface or be limited to a geometrically thin shocked re- 
gion. Thr ough ti me differences of the optical/near-UV 
ingresses, l.ll\ 10 estimated that the stand-off distance 
between the shock and the center of the planet is about 
4.2 R p , leading to a Mg density riM S ^ 400 cm -3 . Here, 
we convert this minimum density to hydrogen number 
density by u sing t he observed metallicity of WASP-12, 
[M/H] = 0.3 (|H09h 



= io( e M E -«H) 10 o.3 



6.76 x 10" 



(2) 



where 6M g = 7.53 and e# = 12.00 are th e solar element 
abun dances for Mg and H, respectively (|Grevesse et al.l 
12001 . Therefore, the minimum density of the shocked 
stellar coronal plasma is 



400 cm" 3 
6.76 x 10- £ 



6 x 10 6 cm" 3 . 



(3) 



The density behind such a shock would, in the adiabatic 
limit, be at most four times the density ahead of the 
shock. From Equation ([3]), this would then require a 
pre-shock coronal density of 



^obs 



1.5 x 10 6 cm -3 



(4) 



in the external ambient medium at the orbital radius of 
the planet. 

The temperature of the absorbing material must be 
such as to allow for the presence of Mg II lines. For 
an adiabatic shock, the temperature immediately behind 
the shock depends on the squared Mach number of the 
flow impacting on the planet. To determine the temper- 
ature and ionization structures of the atmosphere of the 
planet, and therefore quantitatively evaluate the optical 
depth of Mg II lines, detailed radiative transfer models 
are needed. 

3. HYDROSTATIC CORONA 

In this section, we assume that the orbit of the planet 
lies within the stellar magnetosphere. In this case, the 
planet moves through a medium that is confined by a 
rigid stellar magnetic field and so rotates at the stellar 
rotation rate. The relative azimuthal velocity between 
the planet and this medium is therefore 



Am = \uk 



where u„ 



4 ^,cor | 



R, 



orb 



1/2 



2vri? 



P* 



(5) 



= 27ri? or b/-P* is the velocity of the medium 
corotating with the star, and P* is the stellar period of 
rotation. If Au > c s . a bow shock forms around the 



Early UV Ingress in WASP-12b: Measuring Planetary Magnetic Fields 



3 



(a) 

dayside-shock 
(6=90°) 



(b) 

ahead-shock 
(9=0°) 



Au-0 



u,~0 



Aulln 



(c) 

intermediate case 
(6=atan(u r /Au)) 



\ 1 / 



Au 



Fig. 1. — Sketch of shock types (not to scale): (a) dayside-shock (9 = 90°), (b) ahead-shock (9 = 0°), and (c) intermediate case. Arrows 
radially leaving the star depict the stellar wind, dashed semi-circles represent the orbital path, 9 is the deflection angle between n = Au — u r 
and the relative azimuthal velocity of the planet Au. 

planet. We may therefore define a critical stellar rotation 
period such that Au = c s and 



P-, 



2irR, 



orb 



\UK 



(G) 



Figure [2] presents P cr it,* as a function of the coronal 
temperature T. The shaded areas show regions on the 
Pcrit,* — T parameter space where no shock is formed. In 
our notation, P cr it,* < refers to cases where the star is 
counterrotating with re spect to the planet ary orbital ro- 
tation (e.g., WASP-8b: iQ.ueloz et al.ll2010fl . The vertical 
line in Figure [5] represents a critical coronal temperature 
where uk = c s and |P C rit *| - ► °o (cf. Equation (J6))). For 
WASP-12, we obtain 



GJVL 



= •>■■■). I :l0 6 <yA f.^ M ^ t ! K = 4.16x10 s K, 



(Porb/P©) 



(7) 

where m = fim p is the mean particle mass, ks is the 
Boltzmann constant, and m p is the proton mass. We 
adopt fi — 0.66. Above this temperature, no shock for- 
mation is possible at any positive P*. For all the re- 
maining areas of the plot an ahead-shock will develop. 1 
Althoug h P, is unknown, W ASP-12 is likely to be a slow 
rotator (jFossati et a l. 2010b). With this information, the 
presence of an ahead-shock for WASP- 12b constrains the 
coronal temperature to be T < (4 — 5) x 10 6 K. 
For an isothermal corona, the density is 



n(R 



uq exp 



GM*/R* ( R 



ksT/m 
k B T/m \ P 2 



Porb 
p2 



1 



(8) 



where no is the density at the coronal base. From Equa- 
tion ((U), we note that the more slowly the star rotates, 
the lower is the density at any given radius. For the 
very slow stellar rotation rates expected for stars with 
detected transiting planets, P* has little influence on the 
density. The solid lines in Figure |3] show the density 
variation for a range of orbital radius for three differ- 
ent temperatures and a solar value P* = 26 d. The 

1 A "behind-shock" is formed when uk — Wip.cor < 0, or 
Petit * < Porb — Id for WASP-12b. In such situation, the planet 
orbits beyond the Keplerian corotation radius. The coronal plasma 
lags behind the planetary motion and the shock trails the planet. 
Observationally, we would detect a late egress, instead of an early 
ingress. However, Figure ^ shows that for P C rit * < 1 d (WASP- 
12b), the relative azimuthal velocity Au is subsonic, therefore, no 
shock will be formed. 
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Fig. 2. — Critical stellar rotation period (Equation (O) required 
for the formation of a bow shock as a function of the environment 
temperature. Negative periods refer to cases where the planet is in 
a retrograde orbit. 

density was scaled to match the observationally derived 
value of n bs — 1-5 x 10 6 cm~ 3 (Equation (|4])), which is 
the minimum density required for the detection of the 
early ingress. Except for the very low temperature case 
(black solid line), a typical solar coronal base density 
(i.e., uq ~ 10 8 cm -3 ; Withbroe 1988) provides the con- 
dition for the existence of a shock. 

We note that the observationally derived stand-off dis- 
tance from the shock to the center of the planet can 
be taken as approximately t he exten t of the planetary 
magnetosphere tm ~ 4.2 R p ()LHV10l ). Pressure balance 
between the coronal total pressure and the planet total 
pressure requires that, at tm, 



p c Au 2 B c {R mh f 



47T 



B p (r M y 

47T 



■Pp, 



(9) 



where p c — mn b s , p c and P c (P rb) are the local coronal 
mass density, thermal pressure, and magnetic field inten- 
sity, and p p and B p (tm) are the planet thermal pressure 
and magnetic field intensity at r^j- Neglecting the ki- 
netic term and the thermal pressures 2 , we may therefore 

2 It is straightforward, however, to show that the kinetic term 
and both thermal pressure terms are negligible relative to the mag- 
netic pressure terms. From the estimated planetary magnetic field 
(24 G), the magnetic pressure is pg ~ i? p (rjv/) 2 /(47r) ~ 8 X 
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Fig. 3. — Coronal density of an isothermal plasma for three dif- 
ferent temperatures. 



rewrite previous equation as 



B c {Rorb) — B p (r M ), 



(10) 



from which we can estimate an upper limit for the plan- 
etary magnetic field intensity. 
For a stellar dipolar magnetic field, B c (R OT b) = 

-B* (i?*/i? rb) 3 - Assuming that the planetary magnetic 

field is also dipolar, then B p (tm) — B p (R p /rM) 3 ■ P* 
and Bp are the magnetic field intensities at the stellar and 
planetary surfaces, respectively. From Equation (jlOl) . we 
have 



B n — P* 



R*/Rq 

R P /r 



— P* 



1/3.15 
1/4.2 



2.4P*, (11) 



where we have used the observationally derived char- 
acteristics for the WASP-12 system: ru — 4.2 R p 
and i? or b = 3.15 P*. So far, tentative measure- 
ments of the stellar magnetic field have not pro- 
vid ed significant detections (for theoretical p redictions, 
see [Christcnse n. Holzwarth. fc Reineri [2009) , but have 
yielded an upper limit of 10 G in the longitudinal com- 
ponent (|Fossati et al.ll2010H) . Adopting this upper limit 
as P», our model predicts a maximum planetary mag- 
netic field of about 24 G. 

4. STELLAR WIND 

In this section, we investigate the scenario where the 
stellar magnetic field is not strong enough to confine 
the coronal plasma, which expands in the form of stellar 
wind. For a hot corona , the therma l pressure gradient 
is able to drive a wind (|Parker|[l958f) . We consider the 

10~ 3 dyn cm 2 . The kinetic term p c Au 2 /2 in Equation J9j is de- 
pendent on the period of the star, which for WASP-12 is unknown. 
For example, for P* = 26 d, p c A« 2 /2 ~ 4 X 10 -4 dyn cm -2 . 
For T = 2 X 10 6 K, p c ~ 4 X 10 -4 dyn cm -2 . To estimate 
the thermal pressure of the planet, we use values derived from 
Murray-Clay ct al. (2009) for planetary density (~ 10 7 cm -3 ) and 
temperature (10 3 ' 5 K), which result in p p ~ 4 X 10 -6 dyn cm -2 . 
Because pg is much larger than the aforementioned pressures, the 
magnetic terms in Equation J9} will dominate. 



Fig. 4. — Coronal wind velocity of an isothermal wind for three 
different temperatures. Filled circles represent the location of the 
sonic point (u r = c s ) and dashed lines represent the regions of 
supersonic velocities, which, for T = 1 X 10 6 K, is outside the 
range of orbital radii shown above. 



simplest case of a thermally-driven wind which is not 
magnetically channeled. The wind radial velocity u r is 
derived from the integration of the differential equation 



pu r 



du r 
dr 



dp 
dr 



GM* 



(12) 



Figure U presents the wind velocity profile at different 
orbital radius up to a few stellar radii for three isothermal 
wind temperatures. In the distance range P or b < 4P„, 
the wind is still accelerating, and the radial velocities 
are considerably smaller than the terminal velocities 
achieved: Uoo ^ 440, 670 and 1020 km s" 1 for T = 1, 
2, and 4 x 10 6 K, respectively. The filled circles are the 
sound speed of the wind: c s ~ 158 and 223 km s _1 for 
T = 2 and 4 x 10 6 K. Dashed lines represent the regions 
of supersonic velocities. For the lowest-temperature case 
(1 x 10 6 K), the wind becomes supersonic at a much 
larger radii, not shown in Figure |H The vertical dotted 
line shown in Figure |4] represents the orbital radius of 
WASP-12b. 

The density structure of the wind is derived from the 
conservation of mass, where nu r r 2 is a constant of the 
steady-state wind. The dashed lines in Figure [3] present 
the density profile for an isothermal wind, for three coro- 
nal temperatures. In contrast to the hydrostatic corona 
(solid lines), n does not depend on P*. Except for the 
lowest temperature adopted, where the density is well 
described by a hydrostatic density profile, the coronal 
base density no required for the formation of a shock 
in the wind case is larger than the one required by the 
hydrostatic corona case. Nevertheless, as in the static 
scenario, a typical solar coronal density still allows for 
the formation of a shock. 

In the stellar wind case, there are no azimuthal 
forces acting on the flow. Therefore, through conser- 
vation of angular momentum of the particles leaving the 
star (27ri?^/P*), the azimuthal velocity of the wind is 
unwind = 27ri?^/(P*i? or b). This implies that Am for the 
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Fig. 5. — Angle that the shock normal n makes to the relative 
azimuthal velocity of the planet. 



wind case is 
Am = \u K - Up, wind | = 



/ GM* 

V Porb 



1/2 



2nRi 



P*P, 



orb 



(13) 



which is independent of the wind temperature or radial 
velocity. For slow rotators, Au becomes almost indepen- 
dent of P* and is given by Au s» uk- 

The condition for shock formation requires that (u 2 + 
Au 2 ) 1 / 2 > c s , implying that 



IP 



crit,* | 



2rr R 2 JR 



orb 



(u\ + U 2 -C 2 ) 1 /2' 



(14) 



which, qualitatively, produces a similar result as the one 
presented for the static case (Figure [2]). This means 
that the critical temperature derived in the previous 
section (Equation 0) is enhanced by a temperature 
ST = u 2 r m/k B ~ 8 x 10 5 (u r /(100 km s -1 )) 2 K, implying 
that for T < T cr it + ST, a shock will be formed. 

In the wind case, Figure |4] shows that u r ^ at the 
orbital radius of WASP-12b, implying that the shock will 
not form directly ahead of the planet, but at an interme- 
diate angle 8. Because 8 depends on u r (Equation (JT}), 
which depends on the temperature, the orientation of the 
shock is temperature-dependent. Figure [5] shows how the 



angle 9 depends on the wind temperature for P* = 26 d. 
While in the hydrostatic case, we found that an ahead- 
shock (8 = 0°) is always formed, the wind case requires 
a very low temperature to form an ahead-shock. For 
higher wind temperatures, the shock normal gets closer 
to the line of centers between the planet and the star and 
8 approaches 90°. 



5. CONCLUSION 

Motivated by the recent observation on the light 
curve asymmetry of the transit hot- Jupiter WASP- 12b 
(|Fossati et al.ll2010af) , we proposed a model where the in- 
teraction of the stellar plasma with the planet results in 
the formation of a bow shock around the planet, which 
could explain the early transit ingress observed in the 
near-UV. Although we know the orbital radius of the 
planet, we do not know if at this radius the stellar mag- 
netic field is still capable of confining the hot gas of its 
corona, or if this gas is escaping in a wind. Therefore, 
we investigated the physical conditions of the external 
ambient medium around the planet that could allow for 
the formation of such a shock. 

In either case, for plasma temperatures T < (4 — 5) x 
10 6 K we expect that a shock capable of compressing the 
plasma to the observed densities will form around the 
planet. In the case where the coronal plasma is confined, 
we showed that the geometry of this shock is independent 
of the plasma temperature. The shock forms ahead of the 
planet in its orbital path. In the unconfmed (wind) case, 
the angle between the shock normal and the direction of 
planetary motion depends on the ratio of the radial wind 
speed and the azimuthal speed of the planet relative to 
the stellar wind. Since the wind speed depends on the 
temperature, the orientation of the shock is temperature- 
dependent. 

If the planet's orbit takes it through both regions of 
confined plasma and also wind plasma, the orientation 
and density of the shock may change, producing time- 
dependent absorption and duration of transit. 

Given the large range of stellar rotation rates and coro- 
nal/wind temperatures at which a sh ock is c apable of 
producing the densities estimated by ILHV10I . we con- 
clude that this is likely to be a common feature of tran- 
siting systems and may prove to be a useful tool in setting 
limits on planetary magnetic field strengths. For the case 
of WASP-12b, we derived an upper limit for the planet's 
magnetic field of about 24 G. 
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